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Abstract:  Plasmonics  combines  attractive  features  of  nanoelectronics  and  optics  enabling  highly 
integrated,  subwavelength  optical  and  electronic  circuits.  Conventional  metal-based  plasmonic 
devices  have  plasmonic  resonances  less  than  1  pm,  and  thus  metals  interact  effectively  only  with 
ultraviolet  and  visible  light.  The  wide  application  of  plasmonic  devices  hinges  on  practical 
demonstrations  with  low  losses  at  optical  and  infrared  wavelengths  and  this  field  has  received 
much  attention  in  recent  years.  This  review  introduces  different  materials  available  for  plasmonic 
waveguiding  applications  with  an  emphasis  on  highly  conducting  ZnO.  In  addition,  the  paper 
contains  analysis  of  a  set  of  thin  Al-doped  ZnO  (AZO)  layers  grown  by  RF  sputtering  on  quartz 
glass  that  employ  a  unique,  20-nm-thick,  ZnON  buffer  layer,  which  minimizes  the  strong 
thickness  dependence  of  mobility  (p)  on  thickness  (d).  The  values  of  mobility  and  carrier 
concentration  (n)  were  used  to  predict  optical  properties  through  the  Drude  dielectric  function. 
The  optical  and  electrical  properties  of  the  AZO  are  also  used  to  design  insulator-conductor- 
insulator  waveguides  for  long  range  plasmons  using  full-wave  electromagnetic  models  built  with 
finite  element  method  simulations.  In  many  cases,  the  waveguide  has  subwavelength  dimensions 
showing  that  IR  light  can  be  manipulated  in  semiconductor  materials  at  dimensions  below  the 
diffraction  limit. 

Keywords:  microcavity;  polariton;  strong  coupling;  ZnO 
INTRODUCTION 

The  field  of  plasmonics  has  received  much  interest  in  recent  years  and  several  research  efforts 
have  been  aimed  at  the  understanding  of  the  underlying  physics1.  Recent  work  has  focused  on 
demonstrations  that  overcome  the  fabrication  and  material  constraints  and  enable  maturation  of 
the  technology  into  field-testable  applications"  like  optical  computing  and  chips,  enhanced  signal 
detectors,  etc3.  Surface  plasmon  polaritons  (SPPs)  are  quasi-particles  or  excitations  that  result 
from  resonant  coupling  of  photons  to  the  collective  oscillations  of  conduction  electrons  in 
materials4.  SPPs  can  be  described  as  two  dimensional  (2D)  bound  electromagnetic  waves  that 
propagate  along  conductor-dielectric  boundaries5  and  exponentially  decay  away  from  this 
interface5’7.  These  waves  are  characterized  by  strongly  enhanced  localized  fields  and  high  spatial 
confinement,  sometimes  at  subwavelength  scales,  which  makes  them  attractive  for  the  design  of 
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components  not  constrained  by  the  diffraction  limit  which  may  reduce  size,  weight  and  power 
consumption  when  compared  to  conventional  systems1-6. 

While  the  physics  behind  SPPs  was  initially  explored  decades  ago,  the  field  has  experienced 
rapid  growth  in  recent  years  enabled  by  advances  in  fabrication  techniques  which  allow  for 
spatial  resolution  of  nanoscale  features.  Additionally,  the  refinement  of  semiconductor 
processing  and  exploration  of  new  materials,  such  as  highly  conductive  semiconductors9,  has 
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helped  alleviate  losses  that  have  hindered  the  widespread  application  of  conventional  metal- 
based  plasmonic  devices  in  the  optical  and  infrared  (IR)  regime.  Recent  papers  have  shown  that 
alternative  plasmonic  materials  such  as  semiconductor  oxides10'12,  ceramics13,  semimetals14, 
organics15  and  transition-metal  nitrides16  can  offer  advantages  such  as  lower  losses  compared  to 
metals,  and  tunability  of  optical  properties  for  different  regions  of  the  spectrum  including  the 
IR17,18. 

In  this  review,  we  focus  on  transparent  conducting  oxides  (TCOs)  as  potential  candidates  for 
plasmonic  applications  in  the  IR  regime.  Specifically,  we  use  the  example  of  highly  doped  ZnO 
to  demonstrate  the  deposition  methods,  characterization  and  simulation  of  plasmonic  waveguides 
using  experimentally  measured  properties  to  elucidate  the  efficacy  of  these  materials  in  designs 
where  noble  metals  are  the  nonn.  Application  of  metallic  plasmonic  devices  is  relegated  to  the 
ultraviolet  (UV)  and  visible  (VIS)  spectral  ranges  (<1  pm)  since  the  carrier  concentration  (n) 
which  determines  the  plasma/resonance  frequency  (cop  oc  n12)  is  difficult  to  change  in  metals19  20. 
Metals  also  exhibit  large  losses  in  the  IR  region21.  Both  these  issues  can  be  overcome  using 
TCOs  e.g.  ZnO  heavily  doped  with  A1  such  that  n>10  cm'  translating  to  a  resonant  wavelength 
of  approximately  1pm  with  much  lower  losses  than  metals10,22. 

PROGRESS  HISTORY 

Efforts  to  build  plasmonic  devices  that  combine  the  high  bandwidth  of  photonics  with  the 
compactness  of  electronics  ’  "  and  can  be  broadly  divided  into  two  categories  of  interest:  (i) 
Sensing  applications  where  local  enhancement  of  electromagnetic  (EM)  fields  using  careful 
design  of  geometric  properties  of  periodic  structures  such  as  metal-hole  arrays26,27  and  metal 
inclusions"  enables  signal-to-noise  ratio  gains  as  well  as  strong  nonlinear  optical  effects  ’  ;  and 
(ii)  Light  guiding  applications  where  suitable  material  design  and  optimized  guide  geometries 
can  enable  confinement  of  optical  modes  in  subwavelength  waveguides30'32.  This  review 
focusses  on  waveguides  for  long-range  plasmon  excitation  and  propagation.  Traditionally 
waveguides  are  made  of  dielectric  media  and  operate  on  the  principle  of  light  propagation  using 
total  internal  reflection.  The  minimum  dimension  requirement  to  confine  a  guided  photonic 
mode  is  set  by  the  diffraction  limit  and  is  typically  on  the  order  of  Ao/nw  where  Xo  is  the 
wavelength  in  free  space  and  nw  is  the  refractive  index  of  the  waveguide.  In  contrast,  plasmonic 
waveguides  can  sustain  guided  modes  in  subwavelength  waveguides  using  SPPs  at  conductor- 
dielectric  interfaces23,33'35. 

36 

Many  configurations  of  SPP  waveguides  have  been  proposed  to  date  that  leverage 
characteristics  such  as  mode  confinement,  and  low  propagation  or  coupling  losses  .  They  can  be 
broadly  classified  as  metal-insulator-metal  (MIM)  or  insulator-metal-insulator  (IMI) 
waveguides38'40.  The  mode  consists  of  the  SPP  waves  that  travel  along  each  of  the  two  dielectric- 
metal  interfaces.  A  few  examples  of  plasmonic  waveguides,  comparable  to  strip  waveguides  in 
photonics,  are  thin  metal  films41'4  ,  cylindrical  nanowires23,4  '49,  nanostrips  on  dielectric 
substrates36,38,50,51,  dielectric-loaded  SPP  waveguides52'56  and  chains  of  metal  nanoparticles57'60. 
Plasmonic  waveguides  comparable  to  slot  waveguides  are  nanogaps  in  thin  metal  films24,25,39,61, 
channel  plasmon  polariton  waveguides62'65  and  nanogrooves  cut  in  metal  media64'67  e.g.  V- 
grooves  ,  rectangular  nanogaps  in  metallic  surrounding  films  '  ,  nanoholes  in  metallic 
films25,44  and  inverse  slots  formed  by  sharp  metal  wedges42,64,72,73.  While  theoretically  plasmonic 
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waveguide  are  attractive,  each  design  must  be  considered  for  its  individual  strengths  and 
weaknesses  in  the  context  of  fabrication  constraints  and  practical  imperfections  which  inevitably 
diminish  predicted  perfonnance.  For  example,  nanoholes  and  nanoparticle  chains74  are  difficult 
to  fabricate  and  exhibit  dissipative  losses-  .  Strip  waveguides  ’  are  easier  to  fabricate  but 
susceptible  to  scattering  losses  due  to  surface  roughness  and  geometrical  aberrations  .  We 
consider  the  hybrid  configuration,  which  aims  to  combine  the  tight  subdiffraction  confinement 
offered  by  plasmonic  waveguides  with  longer  propagation  lengths  of  dielectric  structure77,78. 
This  can  be  constructed  by  either  embedding  a  conductor  waveguide  in  a  dielectric  media  (IMI) 

70 

or  fonning  dielectric  nanowires  coupled  to  a  metal  (MIM)  .  We  focus  on  the  IMI  hybrid 
waveguide  structure. 

The  presence  of  metal  in  plasmonic  waveguides  introduces  large  ohmic  losses  in  the  IR  regime 
and  has  been  long  considered  a  major  hurdle  to  wide  application  of  these  devices  .  However,  the 
presence  of  a  conducting  medium  in  the  same  spatial  location  as  the  propagating  electric  field 
enables  control  of  mode  characteristics"  ’  .  In  addition,  the  ability  to  exercise  electrical  control 
in  the  same  waveguide  that  carries  optical  information  opens  the  door  to  many  interesting  hybrid 
electro-optical  integrated  components  .  Therefore,  the  best  materials  for  plasmonic  components 
would  minimize  or  eliminate  ohmic  losses  while  maintaining  the  conductivity  that  allows  for  the 
electrical  characteristics19,35.  TCOs  may  be  a  good  alternative  material  to  implement  practical 
devices  that  combine  these  two  features. 

TCOs  have  been  extensively  researched  for  electrodes  and  find  applications  in  liquid  crystal 
displays,  solar  cells  and  light-emitting  diodes  ".  While  the  most  commonly  utilized  TCO  material 
is  indium  tin  oxide  (ITO),  its  high  cost  can  be  prohibitive  and  a  strong  candidate  to  replace  ITO 
is  ZnO  doped  with  Al,  Ga,  or  In.  Al-doped  ZnO  (AZO)  has  also  recently  been  proposed  ~  as  a 
plasmonic  material  in  the  IR  region,  a  range  not  accessible  to  metal-based  plasmonics10,84.  AZO 
can  be  easily  fabricated  into  waveguides  using  standard  fabrication  techniques  e.g. 
photolithography  and  patterned  into  periodic  structures  such  as  conductor-hole  arrays  and 
nanowire  arrays85’87.  We  believe  that  the  material  growth  techniques  and  some  of  the  simulations 
presented  for  waveguiding  applications  can  be  easily  extended  to  periodic  arrays  for  radiation 
enhancement. 

FABRICATION  AND  METHODS 

ZnO  is  an  attractive  material  for  plasmonic  applications  because  high-quality  films  can  be  grown 
on  a  variety  of  substrates  using  many  different  methods  such  as  sputtering  and  pulsed  laser 
deposition  (PLD)89,90.  Material  properties  can  be  tuned  by  controlling  carrier  concentration  (n) 
and,  to  a  lesser  extent,  mobility  (p).  We  have  recently  grown  ZnO  using  PLD  and  simple  air 
anneals  to  red-shift  the  plasmonic  resonance  to  IR  wavelengths10,91.  However,  TCO  films  can 
have  non-unifonnities  related  to  growth  processes  and  the  non-stoichiometric  nature  of  the 
materials  which  causes  optical  properties  to  exhibit  thickness  dependence92’95.  Further,  the 
substrate  interface  can  have  trap  states  which  reduces  overall  carrier  concentration.  These  losses 
are  thickness  dependent  due  to  scattering  observed  at  the  microstructure  level  of  the  films  ’  . 
Another  problem  with  films  grown  on  lattice  mismatched  substrates  is  the  dependence  of 
electrical  properties  including  resistivity  (p),  p,  and  n,  on  the  thickness  of  the  film93,95'97.  It  is 
hypothesized  that  poor  crystallinity  near  the  substrate/film  interface  gives  rise  to  this  problem. 
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Since  unifonn  thin  films  with  low  scattering  losses  are  key  to  producing  practical  plasmonic 
devices,  it  is  not  only  necessary  to  understand  the  origins  of  thickness  dependence  of  optical 
properties  but  also  to  overcome  this  issue  using  careful  design  of  growth  processes  .  Recently, 
Itagaki  et  al"  demonstrated  improved  crystallinity  of  RF-sputtered  ZnO  by  inserting  a  thin 
ZnON  buffer  layer  between  the  substrate  and  ZnO  layer.  For  example,  in  undoped  ZnO  grown 
on  c-plane  AI2O3,  the  rocking-curve  FWHM  of  the  (002)  diffraction  dropped  from  0.490°  to 
0.061°,  and  in  AZO  grown  on  quartz-glass  substrates,  the  thickness  dependences  of  p,  p,  and  n 
was  greatly  minimized.  We  also  observed  low  surface  roughness  and  few  defects  in  these  films 
which  minimizes  losses  . 

We  present  results  of  a  set  of  thin  (d=25-147  nm)  AZO  layers  grown  using  RF  magnetron 
sputtering  at  300°C  in  4/20.5  [seem]  ISb/Ar  ambient,  producing  a  total  pressure  of  0.35  Pa.  For 
each  thickness,  an  AZO  film  was  grown  on  a  bare  quartz  substrate  while  another  was  grown 
under  similar  conditions  on  a  substrate  with  a  20-nm-thick  ZnON  buffer  layer.  The  targets  for  the 
buffers  were  99.99%  pure  ZnO  and  the  applied  RF  power  was  100W.  The  AZO  films  were 
grown  at  200°C  in  pure  Ar  ambient  using  ZnO  targets  with  2wt%AE03.  The  buffer  layers  were 
20-nm  thick,  but  results  not  presented  here,  have  shown  negligible  differences  in  results  for 
buffer  thicknesses  from  5-100  nm.  Electrically,  the  buffer  layers  are  semi-insulating  with 
minimal  effect  on  the  conductivity  of  the  AZO  film.  The  purpose  of  the  N  in  the  buffer  layer  is  to 
inhibit  the  strong  nucleation  tendency  of  ZnO,  which  leads  to  small  grain  sizes  and  results  in 
larger  scattering.  The  buffer  increased  grain  sizes  from  38nm  to  68nm  measured  using  X-ray 
diffraction  co-29  scans.  The  buffer  also  reduces  the  strong  thickness  dependence  of  p  and  n100'102. 

CHARACTERIZATION  OF  ZnO  AS  PLASMONIC  MATERIAL 

Detailed  Hall-Effect  and  reflectance  measurements  were  conducted  to  quantify  the  effect  of  the 
buffer  layer  on  the  quality  of  the  AZO  films  (Figure  1).  Hall-effect  measurements  were  carried 
out  in  a  LakeShore7507  apparatus.  The  magnetic  field  strength  was  IT  and  the  temperature 
range,  6-320  K.  The  samples  were  of  size  1cm  x  1cm,  and  ohmic  indium  dots  were  placed  at  the 
four  corners.  Spectroscopic  reflectance  measurements  were  performed  in  a  Perkin-Elmer900 
UV/Vis  spectrometer  at  294K  over  a  wavelength  range  190-3200  nm.  The  volume  carrier 
concentration  was  independent  of  thickness  and  not  affected  by  the  buffer.  However,  p  was 
strongly  influenced  by  the  buffer,  exhibiting  decreased  thickness  dependence  and  increased 
magnitude  at  a  given  thickness.  The  measured  data  shows  that  n  does  not  depend  on  thickness 
for  either  buffered  or  unbuffered  samples  nor  is  it  affected  by  the  presence  of  the  buffer. 
Additionally,  p  does  depend  on  thickness  and  is  affected  by  the  presence  of  a  buffer.  This  can  be 
explained  by  the  occurrence  of  interface  scattering  where  electrons  in  thinner  samples  scatter 
more  from  the  interface  because  they  are  closer  to  it  on  the  average.  The  buffer  reduces  this 
scattering,  and  also  reduces  or  eliminates  the  trapping  of  free  electrons  from  the  bulk  at  the 
interface97. 

The  measured  electrical  properties  are  used  to  calculate  the  optical  properties  of  the  samples. 

20  3 

High  carrier  concentrations  (n>10“  cm"  )  result  in  large  plasma  frequency  which  in  turn  results  in 
optical  properties  akin  to  Drude  metals  ’  ’  "  .  Samples  of  degenerately  doped 

semiconductors  like  the  samples  presented  here  can  be  described  using  the  Drude  function  that 
describes  pennittivity  as  a  function  of  frequency  (Equation  1).  The  square  root  of  n  directly 
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relates  to  the  plasma  frequency,  and  thus  the  resonance  frequency  can  be  easily  modified  by 
changing  the  carrier  concentration. 


where  sr  and  8;  are  the  real  and  imaginary  part  of  the  permittivity  respectively,  Sint  is  the  intrinsic 
pennittivity  of  the  undoped  material,  nopt  is  the  optical  doping  density  or  carrier  concentration,  co 
is  the  operating  frequency,  e  is  the  electron  charge,  m*  is  the  effective  mass,  and  y  is  the 
relaxation  rate  which  can  be  approximated  by  y=e/poptm*  where  popt  is  the  optical  mobility. 
Comparisons  of  Hall-effect  and  reflectance  measurements  in  these  samples  have  shown  that 
nopt~n,  and  popt~P-  The  permittivity  calculated  using  the  Drude  function  is  plotted  vs.  film 
thicknesses  in  Figure  2  over  a  wavelength  range  of  4-  10pm  for  both  the  buffered  and  unbuffered 
films. 

Both  sets  of  films  show  similar  e,  values  except  at  smaller  thicknesses  where  interference 
scattering  comes  into  play.  The  buffered  films  show  constant  s;  even  in  the  thinner  samples  and 
the  buffer  seems  to  mitigate  scattering  issues.  In  the  unbuffered  samples,  8;  falls  in  the  thinner 
samples  and  this  is  more  pronounced  at  longer  wavelengths.  The  value  of  er  varies  significantly 
depending  on  the  presence/absence  of  the  buffer  layer  but  does  not  change  much  with  thickness 
in  the  buffered  samples.  It  increases  in  magnitude  (negative)  as  wavelength  increases  but  remains 
almost  constant  at  that  value  across  samples  of  different  thickness  in  the  buffered  samples.  The 
biggest  variation  within  these  samples  is  seen  at  the  longest  wavelength  (X=10pm)  across 
thicknesses.  Interestingly,  this  is  not  the  same  trend  for  the  unbuffered  samples  where  we 
observe  that  the  sr  is  unchanged  in  the  thinner  samples  irrespective  of  wavelength  but  drastically 
decreases  with  thickness  and  the  slope  of  the  drop  increases  with  wavelength97.  This  can  be 
attributed  to  the  fact  that  the  mobility  in  the  thinner  films  for  the  buffered  samples  is  larger 
compared  to  the  unbuffered  samples  and  related  to  lower  damping,  lower  losses  and  larger 
crossover  frequencies  where  sr<0  .  The  pennittivity  directly  affects  the  analytical  dispersion 
relation  and  thus  the  waveguide  characteristics  such  as  propagation  loss,  mode  confinement,  etc. 
are  dependent  on  these  quantities106’107.  These  values  are  also  used  in  the  EM  simulation  and  thus 
the  electric  field  associated  with  each  mode,  as  expected,  is  affected  by  the  permittivity  of  the 
AZO  and  the  surrounding  dielectric  film. 

DESIGN  OF  THE  PLASMON  WAVEGUIDE 

Next  we  demonstrate  how  the  optical  properties  obtained  from  the  Drude  model  can  be  used  to 
design  a  practical  waveguide  structure  for  excitation  of  long  range  plasmons,  which  in  many 
cases  supports  propagation  at  the  subwavelength  scale.  The  waveguide  structure  demonstrated 
here  is  comparable  to  the  IMI  structure  in  metals  ’  .  The  geometry  consists  of  an  Insulator- 

TCO-Insulator  where  a  straight  TCO  strip  waveguide  is  embedded  in  a  dielectric  medium.  The 
structure  we  analyze  is  a  symmetric  slab  waveguide  structure  (Figure  3)  consisting  of  a  straight 
rectangular  AZO  strip  waveguide  surrounded  by  a  polymer  called  benzocyclobutene  (BCB) 
which  has  a  refractive  index  of  -1.535.  The  AZO  strip  has  the  following  dimensions:  width=6 
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pm  and  thickness  ranging  from  25-147nm.  The  AZO  film  substitutes  the  metal  in  the  typical 
IMI  setup  and  the  BCB  film  provides  the  insulator  layer. 

Further,  the  waveguide  is  a  hybrid  plasmonic-photonic  structure  where  the  electric  field  travels 
along  the  AZO-polymer  interface  as  an  SPP  mode  and  the  evanescent  energy  that  leaks  from  the 
mode  is  then  confined  by  the  dielectric-air  boundaries  within  the  larger  photonic  guide  fonned 
by  the  polymer.  This  resulting  hybrid  mode  has  lower  propagation  losses  than  a  plasmonic 
waveguide  formed  with  conductor-air  interfaces  alone77.  Plasmonic  waveguide  geometry  and 
materials  can  be  optimized  either  to  maximize  propagation  length  or  confinement  to  the 
conductor  guide.  This  trade-off  is  discussed  in  greater  detail  in  a  later  section.  Also  the  operating 
wavelength  is  critical  to  the  design  since  it  not  only  affects  the  electric  field  due  to  the  geometry 
of  the  waveguide  but  also  because  the  material  parameters  (e.g.  permittivity)  have  strong 
wavelength  dependence.  This  gives  the  scientist  two  disparate  and  uncorrelated  design  choices 
that  can  be  adjusted  to  tailor  and  optimize  the  waveguide  to  best  operate  for  the  specific  need. 
The  first  is  the  geometry  of  the  device  and  the  second  option  that  is  available  with  TCOs  (but  not 
metals)  is  material  parameters  by  changing  doping  levels  and  in  turn  carrier  concentration  as  well 
as  mobility. 

CHARACTERIZATION  OF  WAVEGUIDES 


The  waveguide  is  analyzed  for  propagation  and  confinement  properties  using  an  analytical  model 
while  the  electric  field  profile  is  examined  using  finite  element  method  (FEM)  electromagnetic 
simulations.  The  values  of  the  material  constants  for  the  AZO  layer  in  both  sets  of  models  are 
based  on  experimental  measurements  conducted  on  films  that  were  fabricated  by  Itagaki  et  al. 
The  experimental  values  of  mobility  and  carrier  concentration  were  used  in  the  Drude  model  of 
dielectric  function  to  predict  the  pennittivity  of  the  AZO  film  as  a  function  of  wavelength  and 
thickness,  which  is  subsequently  used  to  predict  the  electric  field  characteristics  using  the  FEM 
model  and  propagating  plasmon  mode  features  including  losses  and  confinement  using  the 
analytical  model. 

First  we  present  an  analytic  model  of  the  waveguide  to  detennine  the  SPP  dispersion, 
propagation  and  localization.  We  assume  the  films  are  smooth  and  lossy  with  wavelength 
dependent  material  parameters  derived  from  experimental  measurements.  The  analytical 
dispersion  results  show  mode  splitting  which  results  in  symmetric  and  antisymmetric  field 
distributions38.  The  analytical  solution  to  the  dispersion  equation  for  the  TCO  film  of  thickness, 
d,  centered  at  z=0.  The  dispersion  relations  are  shown  in  Equation  (2)  and  (3). 
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and  the  dielectric  medium  respectively.  Equations  (2)  and  (3)  represent  coupling  modes  between 
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the  SPPs  and  describe  the  symmetry  of  the  electric  field.  Equation  (2)  describes  the  dispersion 
relation  of  the  higher  energy  mode,  L+,  which  is  antisymmetric  with  respect  to  the  z=0  plane. 
Equation  (3),  on  the  other  hand,  describes  the  dispersion  relation  of  the  lower  energy  mode,  L-, 
which  is  symmetric  about  the  z=0  plane  ’  .  The  complex  wave  vectors  are  calculated  using  the 
Eigen  value  solutions  to  the  wave-vector,  (30,  defined  by  the  dispersion  relations  for  wavelengths 
from  >.=4-10  pm.  This  IR  regime  is  of  interest  particularly  to  enhancing  detector  signal-to-noise 
ratio  by  coupling  to  plasmonic  modes.  These  quantities  were  calculated  using  a  2D 
unconstrained  Nelder-Mead  minimization  algorithm  that  used  an  adaptive  simplex  method.  The 
minimization  criteria  used  was  a  tolerance  of  10'10  nm'1.  We  present  simulation  results  using  the 
analytical  as  well  as  EM  models  for  the  higher  energy  mode  L+,  since  the  propagation  losses  for 
this  mode  are  at  least  an  order  of  magnitude  lower  than  that  of  the  lower  energy  mode.  The 
results  for  L-  were  also  calculated  and  have  the  same  trends  as  L+  and  hence  are  not  presented  in 
the  interest  of  brevity. 

The  Eigen  value  solutions  for  Po  were  used  to  calculate  the  real  and  imaginary  part  of  the 
effective  index  for  both  the  buffered  and  unbuffered  samples  (Figure  4).  The  effective  index  can 
be  represented  mathematically  as  neff  =  p/ko.  Po  and  ko  are  propagation  constants  of  the  guided 
wave  and  wave  propagating  in  free  space,  respectively.  The  real  part  of  the  effective  index  is 
related  to  the  phase  delay  introduced  by  the  waveguide  material  per  unit  length  as  compared  to 
propagation  in  vacuum.  This  quantity  also  directly  affects  the  confinement  of  the  mode  in  the 
waveguide66.  As  the  real  part  of  neff  increases,  the  mode  is  more  confined  to  the  AZO  waveguide 
and  thus  experiences  higher  losses.  The  Eigen  value  solutions  for  the  propagation  constant  were 
also  used  to  derive  values  of  propagation  length  and  skin  depth. 

The  optimization  of  plasmonic  waveguides  to  the  specific  application  of  interest  hinges  on  the 
correct  choice  of  geometrical  properties  and  material  parameters  that  balance  the  trade-off 
between  propagation  length  and  skin  depth  which  relates  to  confinement.  With  this  in  mind,  we 
compare  the  performance  of  the  unbuffered  and  buffered  samples  of  equal  thickness  at  different 
wavelengths.  The  propagation  characteristics  are  analyzed  using  the  propagation  length  (plotted 
in  Figure  5(a)).  Propagation  length  of  the  plasmonic  wave  traveling  along  the  conductor- 
dielectric  interface  is  defined  as  the  1/e  decay  length  of  the  electric  field,  as  Lprop  = 
l/[2 7m (/?)].  The  decay  in  the  electric  field  is  caused  by  the  ohmic  losses  of  the  material  which 
result  in  loss  of  energy  through  heat  dissipation7. 

Next  we  calculate  the  skin  depth  of  each  waveguide  structure  as  a  measure  of  the  confinement  of 
the  guide.  The  SPP  electric  field  skin  depth  is  related  to  the  confinement  of  the  plasmon  mode 
and  can  be  defined  as  the  depth  at  which  the  field  falls  to  1/e  in  the  dielectric  medium,  5=l/|kzi| 
where  i=l  and  2  correspond  to  the  AZO  and  the  BOB  respectively.  The  extent  to  which  the 
evanescent  field  leaks  into  the  dielectric  medium,  perpendicular  to  the  conductor-dielectric 
boundary,  reveals  the  confinement  of  the  plasmonic  mode  within  the  conductor  that  makes  up 
the  plasmonic  waveguide7.  The  skin  depth  is  plotted  in  Figure  5(b)  and  the  trend  is  similar  for 
both  buffered  and  unbuffered  samples.  Skin  depth  is  inversely  related  to  confinement  in  the  TOO 
film.  As  skin  depth  in  the  dielectric  decreases,  the  electric  field  is  more  concentrated  in  the  TOO 
film  and  the  mode  becomes  more  confined.  As  expected  the  trend  for  skin  depth  is  the  same  as 
propagation  length.  Increased  confinement  (reduced  skin  depth)  results  in  greater  propagation 
losses  and  lower  propagation  length.  Thus,  a  trade-off  has  to  be  made  depending  on  the 
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application.  Further,  confinement  can  be  increased  at  longer  wavelengths  by 
optimizing/increasing  waveguide  width.  This  is  not  studied  here,  since  the  parameter  of  interest 
was  thickness  and  the  effect  it  has  on  propagation  characteristics. 

ELECTROMAGNETIC  SIMULATIONS  OF  THE  WAVEGUIDE 

The  perfonnance  of  hybrid  plasmonic  waveguides,  such  as  the  one  we  consider  in  this  paper, 
depends  on  both  the  material  properties  of  the  plasmonic  waveguide  and  surrounding  dielectric 
media  as  well  as  their  geometric  dimensions.  We  study  the  effect  of  changing  the  waveguide 
properties  by  using  experimentally  measured  values  of  material  parameters  for  films  with  and 
without  a  buffer  layer  and  changing  the  thickness  of  the  waveguide.  The  Eigen  value  solutions 
for  the  wave  vector,  (30,  obtained  using  the  analytical  model  were  used  to  calculate  the  real  and 
imaginary  part  of  the  effective  index  which  was  then  used  as  initial  guesses  in  the  mode  solver  in 
the  EM  model.  The  model  was  parametrically  simulated  over  a  mid-IR  wavelength  range  of  4- 
10pm  and  repeated  for  different  film  thicknesses  with  associated  mobility  and  carrier 
concentration  measurements  for  both  the  buffered  and  unbuffered  films.  The  waveguide  was 
simulated  using  full  wave  electromagnetic  simulations  with  finite  element  method  (FEM)  based 
software,  COMSOL  Multiphysics. 

The  waveguide  analyzed  in  this  paper  is  a  symmetric  strip  waveguide  is  an  AZO  strip 
surrounded  by  BCB  which  in  turn  is  surrounded  by  air  for  the  purpose  of  simulations.  Since  the 
structure  is  uniform  along  the  length  of  the  device  (both  in  geometry  and  material),  it  can  be 
completely  characterized  using  a  2D  model  of  the  cross-sectional  plane  perpendicular  to  the 
direction  of  propagation.  This  reduces  computational  overhead  and  increases  speed  considerably 
when  compared  to  a  full  3D  model.  The  model  is  primarily  used  to  analyze  the  spatial  extent  of 
the  electric  field  and  profile  of  the  mode  with  respect  to  the  AZO  thin  film  and  the  surrounding 
dielectric.  While  the  mode  profile  does  not  change  with  length,  the  amplitude  of  the  field 
decreases  to  1/e  as  the  propagation  length  is  approached.  We  also  conducted  a  frequency  domain 
study  to  verify  the  propagation  characteristics  and  those  results  were  found  to  match  the 
analytical  model  as  expected  and  are  hence  not  presented.  The  material  characteristics  of  the 
AZO  waveguide  were  defined  using  a  Drude  model  of  the  complex  permittivity  as  a  function  of 
wavelength  as  described  above. 

The  mesh  used  in  the  model  consisted  of  ~64K  points  over  the  entire  geometry  which 
corresponds  to  a  maximum  element  size  of  k/8.  In  the  regions  where  the  field  was  expected  to 
rapidly  change  (e.g.  dielectric-conductor  boundaries)  a  finer  mesh  was  used.  The  computational 
domain  window  size  was  nominally  set  at  150pm  x  120pm.  The  definition  of  boundary 
conditions  can  affect  the  solution  for  electric  field  in  plasmon  waveguides  since  the  propagating 
modes  are  non-radiative  and  closely  bound  to  its  surface.  We  use  scattering  boundaries  on  the 
outer  air  boundaries  which  are  placed  ~50pm  from  the  plasmon  waveguide  at  a  sufficient 
distance  away  such  that  the  plasmonic  modes  are  undisturbed.  FEM  analysis  uses  Maxwell’s 
equations  in  the  usual  form  to  calculate  the  associated  electric  field  mode  profiles106,107. 

These  Eigen  values  derived  in  the  analytical  model  for  each  waveguide  thickness  over  the 
defined  wavelength  range  were  used  as  initial  guesses  in  the  mode  solver  in  the  finite  element 
simulations.  This  helps  the  Eigen  solver  to  converge  faster  and  more  accurately  to  a  solution  for 

8 


Approved  for  public  release;  distribution  unlimited 


the  electric  field  profiles.  The  Eigen  value  solutions  for  the  propagation  constant  were  compared 
to  the  values  obtained  in  the  analytical  solution  and  were  found  to  agree.  The  results  for  electric 
field  distribution  perpendicular  to  the  direction  of  propagation  are  shown  in  Figure  6  for  the 
buffered  samples  calculated  at  7=4  pm  for  thicknesses  ranging  from  25-147nm.  As  the  film 
thickness  increases,  we  notice  greater  mode  confinement  in  the  AZO  guide  and  the  evanescent 
field  in  the  dielectric  is  minimal.  Similar  trends  and  profiles  are  observed  for  the  unbuffered 
samples  and  those  results  are  not  presented  here. 

CONFINEMENT  -  PROPAGATION  LENGTH  TRADE-OFF 

SPPs  have  been  researched  with  respect  to  different  geometries  and  materials  to  guide 
electromagnetic  energy  at  UV/VIS/IR  frequencies.  The  obvious  comparison  at  these  frequencies 
is  to  dielectric  waveguides  in  terms  of  performance109.  Performance  can  be  defined  in  terms  of 
EM  characteristics  (propagation  length,  losses,  etc.)  and  also  mechanical  properties  (size,  weight, 
etc.).  Dielectric  waveguides  confine  light  in  a  waveguide  core  with  a  higher  refractive  index 
than  the  cladding  region  when  the  propagation  angle  is  greater  than  the  critical  angle.  Plasmonic 
waveguides,  on  the  other  hand,  confine  electromagnetic  energy  at  the  conductor-dielectric 
interfaces  depending  on  the  negative  pennittivity  of  the  conductor  region  which  acts  like  the 
waveguide  core.  Surface  localization  in  SPPs  enables  guiding  light  in  structures  with  dimensions 
much  smaller  than  their  dielectric  counterparts,  sometimes  even  subwavelength  defying  the 
diffraction  limit23,44’60,106.  This  is  particularly  useful  where  high  density  of  components2  with  low 
cross-talk  and  interference  can  lead  to  dramatic  improvements  in  system  performance. 

SPP  propagating  modes  are  concentrated  in  the  conductor  core  with  electric  fields  formed  by  the 
coupling  of  SPPs  at  the  two  metal-dielectric  interfaces.  This  mode  coupling  reduces  losses 
resulting  from  the  interaction  with  the  conductor  but  in  turn  also  reduces  confinement  in  both  the 
vertical  and  lateral  dimensions.  Propagation  lengths  of  ~2.5mm  have  been  demonstrated  for 
subwavelength  waveguides  that  have  confinement  in  2D  and  thin  films  of  finite  width  with 
subwavelength  confinement  in  only  one  dimension  have  shown  propagation  lengths  over  several 
millimeters48,50,57.  Tighter  confinement  of  the  mode  reduces  electric  field  cross-section  and 
causes  the  fields  to  be  more  concentrated  to  the  conductor  region  which  in  turn  decreases 
propagation  length  ’  ’  .  Conversely,  allowing  the  energy  to  escape  into  the  dielectric  medium 
can  alleviate  losses  but  reduces  confinement.  Thus  all  plasmonic  waveguide  configurations  must 
balance  the  trade-off  between  confinement  and  propagation  losses  and  be  optimized  in 
accordance  with  the  importance  of  either  in  the  application33,68.  This  trade-off  can  be  quantified 
by  comparing  the  dielectric  skin  depth  and  propagation  length.  Tighter  confinement  to  the 
conductor  waveguide  is  desirable  in  applications  where  the  mode  characteristics  may  be 
modified  by  changing  the  material  characteristics  of  the  guide.  This  is  examined  in  detail  using 
Figure  7  where  the  propagation  length  is  plotted  with  the  real  part  of  the  effective  index  which 
directly  relates  to  confinement  and  is  inversely  proportional  to  dielectric  skin  depth.  As 
confinement  increases,  propagation  length  decreases.  This  contrasting  trend  can  be  balanced 
using  three  features,  namely  width,  height  and  material  of  the  waveguide. 

We  begin  our  discussion  with  the  thickness  of  the  guide.  Not  all  the  waveguide  structures 
explored  for  plasmonic  mode  propagation  can  be  used  for  subwavelength  confinement  nor  are  all 
supported  modes  localized  in  a  structure  with  dimensions  under  the  diffraction  limit.  Conductor 
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strip  waveguides,  like  the  ones  analyzed  in  this  paper,  are  capable  of  either  short  or  long  range 
SPP  propagation.  As  seen  from  the  results,  reducing  film  thickness  reduces  confinement  in  the 
conductor  and  the  propagating  long  range  plasmonic  modes  are  weakly  localized.  These  modes 
are  subject  to  lower  losses  since  only  a  small  portion  of  the  energy  is  carried  in  the  dissipative 
conductor  while  the  rest  propagates  through  the  lower  loss  dielectric  media.  The  longer 
propagation  lengths  are  useful  for  macroscopic  light  transport.  The  hybrid  guide  can  then  be  used 
to  convert  incident  light  from  a  dielectric  mode  (thin  conductor)  that  propagates  for  several 
millimeters  to  an  SPP  mode  (thick  conductor)  at  a  specific  position  along  the  waveguide  to 
enhance  light-matter  interaction  for  applications  such  as  focusing  or  modulation110'111.  The  SPP 
mode  could  then  be  coupled  back  to  a  dielectric  mode  to  transport  light  to  other  applications  or  a 
detector.  Further,  the  diffused  electric  field  of  a  dielectric  mode  can  interact  effectively  with 
ambient  conditions  and  be  used  for  sensing  small  changes  in  refractive  index  such  as  those 
introduced  by  a  target  analyte  for  applications  such  as  chem/bio  sensing1 12-1 14.  On  the  other  hand, 
increasing  the  thickness  of  the  waveguide  increases  the  localization  of  the  mode.  These  modes 
experience  high  losses  and  are  characterized  by  lower  propagation  distances  but  stronger 
confinement  of  the  electric  field.  In  strongly  confined  modes,  SPP  waveguides  can  be  viewed  as 
light  concentrators  that  enhance  light-matter  interaction  ’  .  The  diminished  size  of  these 
structures  makes  them  very  useful  for  highly  integrated  optical  circuits  and  can  even  be  used  as 
subwavelength  optical  devices  in  techniques  where  long  propagation  lengths  are  not  required  to 
maintain  performance. 

While  the  thickness  of  the  waveguide  determines  confinement  in  the  vertical  dimension,  the 
waveguide  width  can  be  tailored  to  alter  the  mode  characteristics  in  the  lateral  dimension75’115. 
For  a  symmetric  IMI  waveguide  constructed  of  metal,  propagation  lengths  of  several  microns 
have  been  reported  for  waveguide  widths  of  -200  nm.  As  expected,  the  propagation  length 
depends  on  the  width  of  the  waveguide  and  previously  published  results  have  shown  that  as  the 
width  of  the  guide  decreases,  the  number  of  supported  modes  also  falls  until  a  certain  width 
below  which  no  modes  are  guided45’75’108,115.  In  contrast,  waveguides  with  micrometer  widths  can 
support  plasmon  modes  that  have  propagation  lengths  of  ~1  cm  at  telecommunication 
frequencies116.  This  same  trend  is  also  observed  in  asymmetric  waveguide  structures  where  the 
conductor  strip  is  fabricated  on  an  insulating  substrate  which  serves  as  the  lower  cladding  and  air 
acts  as  the  upper  cladding.  Although  the  metal-strip  waveguide  can,  in  general,  support  long 
propagation  lengths  for  large  strip  widths,  it  is  not  as  useful  for  the  purpose  of  subwavelength 
confinement. 

Finally,  the  material  that  constitutes  the  waveguide  can  also  be  altered  to  either  maximize 
propagation  length  or  confinement  as  well  as  add  functionality  to  the  device  that  is  not  available 
using  metals.  The  SPP  mode  can  be  confined  in  the  z-direction  when  kz  is  imaginary  and  (s'c  + 
s'd)  <  0  where  s'c  and  s'd  are  the  permittivity  values  of  the  conductor  and  dielectric  respectively. 
Also,  s'c<0  must  be  satisfied  so  that  kx  is  real  and  the  SPP  can  propagate  in  the  x-direction.  If  s'd 
=  1  (air),  then  s'c  =  (s'int  -  s'Dmde)  <  -1  ,  where  s'int  is  the  intrinsic  pennittivity  of  the  undoped 
material  and  s'omde  represents  the  second  and  third  term  of  equation  (1)  detennined  by  the 
material  parameters  of  the  conductor.  Therefore,  s'Dmde  needs  to  be  large  (s'Dmde  =4.72  for  ZnO) 
which  requires  a  high  carrier  concentration  coupled  with  high  mobility.  High  concentration  leads 
to  a  high  refractive  index  which  produces  good  confinement  of  the  SPP  in  the  conductor  with 
minimal  field  leakage  into  the  dielectric.  On  the  other  hand,  high  concentration  results  in  greater 
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losses10'107.  Plasmonic  waveguides  are  usually  made  of  a  noble  metal  such  as  silver  or  gold.  The 
metal  layer,  which  forms  a  very  integral  part  of  plasmonic  waveguides,  is  also  responsible  for 
inherent  ohmic  losses  from  resistive  heating.  Metals  have  very  high  carrier  concentrations  and 
large  plasma  frequencies  which  results  in  high  values  of  imaginary  part  of  pennittivity  that  is 
responsible  for  loss  ’  ’  .  In  addition,  practical  devices  based  on  metals  have  imperfections  like 
surface  roughness  and  grain  boundaries  which  increase  scattering  that  also  contributes  to 
losses16’18’21. 

Losses  can  be  mitigated  either  through  loss  compensation  using  gain  media  or  by  using  lower 
loss  conductors  instead  of  metals20’117'119.  Some  research  groups  have  explored  adding  optical 
gain  using  additional  active  medium  such  as  pumped  dyes  in  the  dielectric  layer  to  compensate 
for  internal,  radiative  and  total  losses  in  the  propagating  plasmonic  waveguide  “  '  .  Another 

approach  to  lower  losses  is  the  use  of  materials  such  as  transparent  conducting  oxides  or 
transition  metal  nitrides  which  offer  a  viable  alternative  in  systems  where  incorporation  of  gain 
media  may  not  be  easily  achieved.  Conducting  materials  such  as  TCOs  can  be  grown  to  have 
appropriate  carrier  concentrations  and  thus  smaller  losses.  These  materials  have  demonstrated 
significantly  lower  losses  when  compared  to  silver,  the  metal  that  exhibits  the  lowest  losses  in 
the  VIS  and  near-IR  ranges.  The  growth  processes  can  also  be  tweaked  to  reduce  losses  even 
more  in  these  alternative  materials.  For  example,  nitrides  can  be  deposited  as  crystalline  layers 
with  minimal  surface  roughness  on  sapphire  which  reduces  scattering.  Another  example  is  the 
use  of  a  buffer  layer  (made  of  ZnON)  used  with  AZO,  as  discussed  in  detail  in  this  paper.  The  N 
in  the  buffer  layer  inhibits  the  strong  nucleation  tendency  of  ZnO,  which  leads  to  small  grain 
sizes.  X-ray  diffraction  co-29  scans  indicate  an  increase  in  grain  size  from  38  to  68  mn  due  to  the 
buffer.  Further  details  on  growth  and  structural  analysis  can  be  found  in  Refs  xx. 

Replacing  metals  in  plasmonic  devices  also  open  the  door  to  additional  functionality  such  as 
static  or  dynamic  tunability.  Conventional  metals  are  not  tunable  and  devices  have  to  be 
designed  to  work  at  a  specific  frequency;  therefore,  multiple  designs  are  needed  for  different 
operating  wavelengths,  which  is  not  ideal  for  multispectral  application.  TCOs  and  other  nitrides 
allow  for  tunability  and  bandgap  engineering  (e.g.,  alloying  ZnO  with  Cd  and  Mg  )  by 
changing  fabrication  processes126’127.  Furthermore,  these  materials  can  be  deposited  using 
techniques  such  as  chemical  vapor  deposition,  atomic  layer  deposition  and  molecular  beam 
epitaxy  which  are  CMOS-compatible  and  may  lead  to  integrated  monolithic  devices. 

CONCLUSIONS 

Recent  research  efforts  have  focused  on  the  expansion  of  structures  as  well  as  material  systems 
that  can  be  applied  in  the  field  of  plasmonics.  These  efforts  have  led  to  the  successful 
demonstration  of  many  plasmonic  structures  that  are  capable  of  guiding  and  manipulating  light 
below  the  diffraction  limit.  We  reviewed  the  progress  of  plasmonic  waveguide  structures  with  an 
emphasis  on  transparent  conducting  oxides  particularly  ZnO  as  an  alternative  material  for  the 
excitation  and  propagation  of  long  range  plasmons.  We  also  discussed  in  detail  the  trade-off 
presented  by  plasmonic  waveguides  between  propagation  length  and  confinement  and  attributes 
of  the  waveguide  such  as  choice  of  materials  and  geometry  that  can  be  used  to  balance  this  trade 
space. 
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Within  this  context,  we  presented  simulation  results  of  Al-doped  ZnO  plasmon  waveguide 
structure  surrounded  by  a  dielectric  medium  that  was  based  on  experimentally  measured  material 
properties.  The  pennittivity,  propagation  loss  and  skin  depth  of  the  waveguide  were  then 
characterized  as  a  function  of  Al-doped  ZnO  film  thickness  (25-147  mn)  and  operating 
wavelength  (4-10  pm)  for  films  grown  both  with  and  without  a  buffer  layer.  The  buffer  layer  was 
shown  to  be  effective  in  reducing  the  thickness  dependence  of  material  properties  of  the  AZO 
layers  and  thus  help  to  improve  propagation  characteristics.  We  achieved  reasonable  values  of 
effective  index  and  a  low  propagation  loss,  which  verifies  that  this  configuration  is  suited  to  the 
excitation  and  propagation  of  long  range  SPPs  in  waveguide  applications.  This  indicates  that 
ZnO  films  grown  at  different  thicknesses  with  buffer  layers  still  exhibit  very  similar 
characteristics  in  plasmonic  waveguide  applications.  This  will  help  overcome  the  discrepancy 
between  simulation  and  experimental  results,  where  fabrication  deviations  in  thickness  can  create 
large  changes  in  quality  of  modes  obtained. 

The  unique  ability  of  plasmonic  components  ~  to  deliver  electromagnetic  energy  via  light  to 
both  optical  and  electronic  devices  makes  them  attractive  for  the  development  of  dense 
integrated  optical  processing  circuits  and  systems.  Since  both  electrical  and  optical  signals  can  be 
transceived  using  the  same  structures,  plasmonic  components  will  likely  provide  the  link 
between  conventional  optical  and  nanoelectronic  components  for  infonnation  and  signal 
processing.  The  key  to  practical  systems  based  on  plasmonic  structures  hinges  on  the 
development  of  theoretical129’130,  fabrication  and  experimental  methods48  that  can  reliably  and 
efficiently,  generate,  manipulate  and  detect  plasmonic  waves  .  Demonstrations  of  simple 

132 

plasmonic  structures  have  inspired  the  recent  development  of  more  complicated  elements  such 
as  resonators  ’  ,  interferometers  ,  splitters  and  integrated  active  devices  that  incorporate 

plasmonics  into  sources,  waveguides  and  detectors.  Plasmonics  holds  the  promise  to 
dramatically  improve  a  wide  number  of  applications  such  as  lab-on-a-chip116’135,  optical  signal 
processing,  sensing  and  imaging136.  We  believe  that  plasmonics  will  enable  optical  and 
electronic  signal  management  in  unprecedented  ways  by  unifying  the  most  desirable  features  of 
electronic  and  optical  devices  in  novel  technology  platforms  with  high  performance. 
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Figure  1 :  Resistivity  (p),  mobility  (p),  and  carrier  concentration  (n)  as  a  function  of  film 
thickness  (d)  for  the:  (a)  Unbuffered  samples;  n  and  p  are  normalized  to  electrical  thickness, 
dei=d-13  mn;  i.e.,  n=ns/(d-13);  p=ps*(d-13)  and  (b)  Buffered  samples;  n  and  p  are  normalized  to 
electrical  thickness,  dei=d+4  mn;  i.e.,  n=ns/(d+4);  p=ps*(d+4).  Reprinted  from  [24]. 
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(b)Permittivity  (imag)  vs.  Thickness 
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Figure  2:  Plot  of  permittivity  (s)  vs.  waveguide  thickness  for  a  symmetric  guide  with  ZnO 
waveguide  embedded  in  BCB.  The  different  color  plots  represent  each  operating  wavelength 
from  A=4-l  0pm.  The  solid  and  dashed  lines  represent  the  buffered  and  unbuffered  samples 

respectively. 
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Figure  3:  (a)  Sketch  of  the  IMI  waveguide  architecture  with  quartz  substrate,  ZnO  waveguide 
and  polymer  BCB  cladding;  (b)  Waveguide  structure  in  cross-sectional  view  with  dimensions  of 

the  guide 


(a)  Effective  index  (real)  vs.  Thickness 


(b)  Effective  index  (imag)  vs.  Thickness 
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Figure  4:  Plot  of  effective  index  (neff)  vs.  waveguide  thickness  for  a  symmetric  guide  with  ZnO 
waveguide  embedded  in  BCB.  The  different  color  plots  represent  each  operating  wavelength 

from  k=4- 1 0pm 
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(a)  Propagation  length  vs.  Thickness 


(b)Skin  depth  vs.  Thickness 
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Figure  5:  Plot  of  (a)  propagation  length  (Lprop)  and  (b)  skin  depth  (5)  vs.  waveguide  thickness 
for  a  symmetric  guide  with  quartz  substrate,  ZnO  waveguide  embedded  in  BCB.  The  different 
color  plots  represent  each  operating  wavelength  from  k=4- 1 0pm 
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Figure  6:  Normalized  transverse  electric  field  intensity  for  symmetric  slab  insulator-TCO- 
insulator  waveguides  at  A,=4pm  for  the  buffered  samples.  Plot  of  electric  field  for  ZnO 
waveguide  embedded  in  BCB  with  ZnO  thickness  of:  (a)  25nm;  (b)  41nm;  (c)  55nm;  (d)  70nm; 

(e)  85nm;  (f)  103nm;  (g)  127nm;  (h)  147nm 
16 

Approved  for  public  release;  distribution  unlimited 


Prop,  length,  Eff.  index(real)  vs.  Thickness 


Figure  7:  Propagation  length  (Lprop)  and  real  part  of  effective  index  (Re(neff))  vs.  waveguide 
thickness  for  k=4pm.  The  electric  field  penetrates  deeper  into  the  BCB  layer  as  the  real  part  of 
effective  mode  index  Re(neff)  approaches  the  refractive  index  of  BCB  (1.535).  As  confinement 

decreases,  propagation  length  increases. 
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